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Abstract 

In this paper we consider a Schrodinger eigenvalue problem with a po- 
tential consisting of a periodic part together with a compactly supported 
defect potential. Such problems arise as models in condensed matter to 
describe color in crystals as well as in engineering to describe optical pho- 
tonic structures. We are interested in studying the existence of point 
eigenvalues in gaps in the essential spectrum, and in particular in count- 
ing the number of such eigenvalues. We use a homotopy argument in 
the width of the potential to count the eigenvalues as they are created. 
As a consequence of this we prove the following significant generalization 
of Zheludev's theorem: the number of point eigenvalues in a gap in the 
essential spectrum is exactly 1 for sufficiently large gap number unless a 
certain Diophantine approximation problem has solutions, in which case 
there exists a subsequence of gaps containing 0, 1, or 2 eigenvalues. We 
state some conditions under which the solvability of the Diophantine ap- 
proximation problem can be established. 

1 Introduction 

Periodic media with defects have long been studied in physics and applied math- 
ematics. Such problems have been studied in quantum mechanics and condensed 
matter physics[T] to model the effect of impurities on the color of a crystal. Since 
the development of photonic crystals [101 [TH [T^] and the ability of engineers to 
construct such defects in an optical context this area has received new attention. 
Most of these studies are concerned with the characterization of the spectrum of 
the medium due to the presence of the defect. More precisely, since a defect does 
not change the essential spectrum of the associated operators, the main problem 
is to study the creation of isolated point eigenvalues of finite multiplicity in the 
spectral gaps of the unperturbed problem. 

There is a fairly extensive literature on the study of the existence and number 
of bound states with energies that lie in a spectral gap. Much of this literature 
[21 [TTl [T^ [TB] employs some variant of the Birman-Schwinger principle in order 
to count the number of gap modes. The Birman-Schwinger principle is, in 
essence, a homotopy argument in the coupling strength of the defect potential, 
where one tries to estimate the number of eigenvalues below a given value of the 
coupling strength. Weyl asymptotics and Dirichlet-Neumann bracketing have 
been used to produce asymptotic estimates for the number of defect eigenvalues 
in the spectral gap of Schrodinger operators in the large coupling limit. Other 
studies fS] have used formal perturbation methods to investigate the origin and 
distribution of defect modes. 

Perhaps the closest results to the ones presented here are due to Zheludev[39l 
liOj . Rofe-Betekovpn 1^ and Firsova[T31 HI] (see also the paper of Gesztesy 
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and Simon [IB], which gives a self contained derivation of these results). In 
these papers the following results are shown: given Schrodinger operator in one 
dimension with potential of the form qper{x) + qdefix), where qperix + a) = 
qperix) is periodic, qper e Ll^iR), (1 + \x\)qdefix) e Li(R) then 

i There are a finite number of eigenvalues in each gap in the essential spec- 
trum. 

ii There are at most two eigenvalues in each sufficiently large numbered gap. 

iii If J qdef{x)dx ^ there is precisely one eigenvalue in each sufficiently 
large numbered gap. 

In this paper we will be concerned with the periodic Schrodinger eigenvalue 
problem with a compactly supported defect 

— Uxx -'r q{x) u ~ E u a; e R (1) 

{qper{x) X < 
qdef{x) a:e(0,l) (2) 
qper{x - 1) a; > 1 

where the potential qp^r has period a : qper(x + a) = qp^rix). Note that we 
do not necessarily assume any particular relationship between the width of the 
defect (scaled for convenience to be 1) and the period of the potential qper{x) 
(denoted by a) so there is a relative phase shift or dislocation in the potential 
between — oo and +oo if is not an integer. In engineering applications 
defects are created by cutting and layering materials with different properties, 
so the above form of the defect is actually extremely natural from the point of 
view of applications [57]. This asymptotic phase shift in the periodic potential 
will turn out to have very interesting physical and mathematical properties, and 
does not appear to have been previously considered. 

An outline of the paper is as follows: we first give a simple topological argu- 
ment relating the number of point eigenvalues in a gap to the number of con- 
nected components of the intersection of the gap with a certain set, essentially 
a resolvent set associated with the defect. Next we give necessary Diophantine 
condition for large- numbered gaps to intersect this resolvent set. In the end 
we establish the following result: Defining Aq to be the mean energy difference 
between the defect and periodic potentials 

= / qdef{x)dx 1 qper{x)dx 
Jo Jo 

there is a set JFq related to the distribution of errors in Diophantine approxi- 
mation of a such that the following holds: if aAq ^ Ta then every sufficiently 
large numbered gap contains exactly one eigenvalue. If aAq € Ta there exists 
a subsequence of gaps related to rational approximants of a which may contain 
0, 1 or 2 eigenvalues. The set Ta satisfies Ta = ^a' if a and a' are related by 
an element of the modular group PSL(2,Z). The results of Zhcludev, Firsova 
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and Rofe-Betekov are analogous to the case with no asymptotic phase shift, 
= mod 1, whose equivalence class is the rationals. For rational a the 
set J-a is simply the point 0, analogous to point (in) above. For irrational a, 
however, the set Fa is much more interesting. We characterize J-a for certain 
cases, and give an explicit example of a periodic potential and defect having a 
subsequence of gaps containing exactly 2 eigenvalues. 



2 Preliminaries 

In this section we define the important quantities for the defect eigenvalue cal- 
culation and present some technical results on the behavior of these quantities. 
Many of these are standard results which can be found in (for example) the text 
of Magnus and Winkler on Hill's equation. 



2.1 Basic Definitions and Equations 

We will usually write equation ([1]) as a system 



d ( u\_ f E-q{x) 



dx \ P 



J 



1 

-1 



where the potential q is assumed to be piecewise . This assumption could 
probably be weakened to allow Hamiltonians defined as quadratic forms, but 
we have not done this. 

It is clear from a standard Weyl sequence argument that the essential spec- 
trum of llj is the same as the essential spectrum of the purely periodic problem 

-Uxx + qper{x)u = Eu 

and, because the defect is compactly supported, there are no imbedded eigenval- 
ues in the essential spectrum. See the general discussion in Reed and Simon [25] 
and more specifically Theorem 6.28 in Rofe-Beketov and Kholkin^^, which ap- 
plies immediately to this case. Thus we need only understand the presence of 
point spectrum in the resolvent set of the periodic problem - the spectral gaps. 
Let us recall some standard definitions and facts from Floquet theory. 

We define M(E), the monodromy matrix for the periodic problem, to be 

M{E) = U{E,a), (4) 

where the fundamental solution matrix U {E, x) satisfies 

U, = j(^^' 'J"^''^'') \ ^U = JHp^riE, x)U (5) 
C/(S,0) = /. 
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where / is the identity matrix. The monodromy M{E) satisfies the identity 
M*JM = J (i.e. M{E) is a symplectic matrix). The Floquet discriminant 
associated to the periodic problem is defined to be 

k{E) = tr(M(£;)) 

The band edges {Ej}J^Q are the roots of k'^{E) — 4 in increasing order and 
respecting multiplicity. The band edges Ek are periodic points {k{Ej) = 2) if 
j = 0,3 mod 4 and are antiperiodic points {k{Ej) = —2) if j = 1,2 mod 4. 
An energy E is said to belong to a band or a gap if k'^{E) < 4 or k'^{E) > 4 
respectively. The bands are given by intervals of the form 

fc2(S)<4 E eiE2,,E2j+i) 

and the gaps by the intervals of the form 

k^iE)>A EeiE2j^i,E2j) = g,. 

The band edges satisfy E2j < E2j+i < £^2^+2, so that the bands are never 
empty, although the gaps may be so. The sign of the derivative of the Floquet 
discriminant k'{E) is fixed on bands: k'{E) ^ for energies in a band, E S 
{E2j , E2j+i), and it can only vanish at a band edge if £'2^+1 = £2^+2 and the 
gap is empty. An empty gap is referred to as a double point, and the monodromy 
matrix M(E) is a multiple of the identity if and only if i? is a double-point. 
The slope of the Floquet discriminant in a band can be interpreted as the Krein 
signature[20','38', T of the eigenvalues of AI (E), although we will not exploit this 
here. The essential spectrum is the union of the bands and band edges: 

Cess = U°^o[-^2j", £'2j+l] 

For energies in a gap, E G (i?2j-i, £'2j), the monodromy matrix M{E) has two 
distinct real eigenvalues, 

M(E)v±=\±v± |A_|<1<|A+| 

one of modulus greater than one and one of modulus less than one. The resolvent 
set of the periodic problem 

will prove important. Here we use the convention that E^i = — cx) so the zeroth 
gap is the semi-infinite one (—00, Eq). 

Similarly N[E, x) will denote the fundamental solution matrix associated to 
the defect 

N{E, 0) = /. 
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We will define the Floquet discriminant for the defect to be 

kdef{E,x) = tr {N{E,x)) 

We will refer to a point {E, x) as belonging to a band, band-edge or gap for the 
defect problem if k'^^j:{E, x) < 4, k'^^^^E, x) = 4, or k1^f{E, x) > A respectively. 
The eigenvectors of N{x, E) in a gap will be denoted 'w± : 

Nw± — T±w± |r_| < 1 < |t+|. 

We will often work with N{E, 1), the fundamental matrix evaluated at the 
end of the support of the defect potential, which with some small abuse of 
notation we refer to as the monodromy map for the defect. This will be denoted 
by N{E) = N{E, 1). We will also have occasion to fix either x or E, and consider 
the the band-gap structure as the other quantity varies. In such cases we will 
refer to the gaps as being "in x" or "in E" as appropriate. 

For the remainder of this discussion we consider a fixed energy E S Gj in a 
gap of the essential spectrum. The eigenvectors v± represent Jost-type solutions: 
v+ represents a solution to the periodic problem satisfying hm„^_oo M"-v+ = 0, 
while represents a solution to the periodic problem which satisfies lim„^+oo Af "tj- = 
0. In order to have an solution the solution which is decaying at —oo must 
connect to the solution which is decaying at +oo. Since the defect potential is 
compactly supported this eigenvalue condition reduces to the following: A'^, the 
transfer matrix for the defect mode, should map v+, the left Jost eigenvector 
for the periodic problem to a scalar multiple of v_, the right Jost eigenvector 
for the periodic problem. In other words for E E Qj C Rper and N{E) defined 
as above, the eigenvalue condition becomes 

N{E)v+^^iv^. (7) 

Since we are in the 2x2 case this is equivalent to the vanishing of the following 
inner product: 

{v_,JN{E)v+) = f{E)^Q, (8) 

where J is the standard skew-symmetric matrix J = ^ ^ ^ defined previ- 
ously. The function f{E) is the Evans function associated to the defect poten- 
tial. Note that, by continuity f(E), can be defined uniquely at the band-edge: 
standard results (see Katof20l) guarantee that the eigenvectors have a Puiseux 
series representation in powers of {E — Ej)^ in the neighborhood of a band edge 
Ej . In a more general setting the Evans function is defined as a determinant of 
exponentially decaying Jost-type solutions, but throughout this paper we will 
make extensive use of the fact (obviously unique to the second order case) that 
the determinant can be written as the inner product with the standard skew- 
form J. This, together with the fact that the Hamiltonian nature of the problem 
respects the skew-form, will simplify many of the calculations. 

The zeros of the Evans function defined in Eq. ([8]) define the defect eigen- 
values. However it is convenient to consider the width of the defect 
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parameter. In this way we can make a homotopy argument in x, and count the 
roots of the above equation as they are created, with a; = corresponding to 
a purely periodic problem and a; = 1 to the periodic problem with a defect. If 
N{E, x) is the fundamental solution operator for the defect potential as defined 
in ^ we define a generalized Evans function to be 



The goal is to apply an implicit function argument to show that the zeros of 
the generalized Evans function lie on continuous curves E{x), and thus can be 
tracked as functions of the homotopy parameter x. 

Throughout this paper Qj — (-^2^+1, £■2^+2) will denote a gap in the essential 
spectrum, which is assumed to be non-empty. We will let Dg. represent the 
number of defect eigenvalues in Qj : 



Of course the goal of the paper is to understand Dg for a general potential. 
Finally TZdef will denote the set 



the energies for which the map across the defect has real eigenvalues. If one 
considers the periodic Schrodinger operator where the potential is the defect 
potential qdef which is defined on [0, 1] is extended to a 1-periodic function in the 
obvious way, then the set TZdef is the union of the resolvent set of this operator 
with the periodic and antiperiodic eigenvalues. Our count of the number of point 
eigenvalues in a gap can be expressed in terms of the number of components of 
the intersection of this set with the resolvent set of the periodic problem. Note 
that, despite the similarity in notation, TZdef is a closed set whereas TZper is an 
open set, for reasons that should become clear. 

2.2 Preliminary lemmas 

We start with some facts about the dependence of the eigenvalues and eigenvec- 
tors of a 2 X 2 real symplectic matrix M{E) on a real parameter E. We assume 
that E £ TZper, so that the eigenvalues are real and distinct. The eigenvectors 
v±, while not orthogonal, are linearly independent and form a basis. We will 
normalize these vectors so that ||w+|| = \\v-\\ — 1. It is convenient to define a 
second set of vectors, the dual basis, J*, J*. It is clear that these two sets 
are bi-orthogonal (though typically not orthonormal). 

Lemma 1. Suppose that M(E) is defined as in ^ has real distinct eigenvalues 
|A_| < 1 < |A+| and eigenvectors v± (\v±\ — 1), and define a± — {^j^, Jv±). 
Then we have 



f[E,x)^{v-,JN{E,x)v+), 



X e [0,1]. 



(9) 



Dg, - #{E e g, I f{E) = 0). 



TZdef{E\\tr{N{l,E))\^>4}, 



a± = 




(10) 
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Similarly if N(E,x) defined in ^ has real distinct eigenvalues t± with |t_| < 
1 < |t+| and associated eigenvectors w± (\w± \ = 1) and we define 

{^,Jw±)^f3^ (11) 

{'^,J^±)=3± (12) 
then the functions (3±,S± satisfy 

P± = ^±<^±^^^y^±) (13) 

d± = -±i-^M^±) (14) 

Note that the normalization \\w±\\ = 1 implies that 

dw± 



dE 
dw± 



P±Jw±, (15) 
5±Jw±, (16) 



dx 

and similarly v± . 

Proof. The eigenvalue equation is 

Mw+ = A+w+, (17) 
Differentiation with respect to the parameter E gives 

dM , dX+ , dv+ .^^ 

Taking the second equation and multiplying by J and dotting with v±, and 
using the fact that M is symplectic, gives 

a+{X^^ - A+) = {v+, J Me v+) 
a_(Al^-A_) = {v^,JMev^) 

Multiplying by J and dotting with Vzp gives the equation for the change in the 
eigenvalue, and similar calculations give the equations for the other parameters. 
Note that all of the denominators are non-zero away from the band edges. □ 

Remark 1. This is a slight modification of the Hellmann-Feynman lemma for 
self-adjoint matrices to the case of matrices in S'p(2,R) with real distinct eigen- 
values. Note that this result only uses the fact that the matrix M is in 5*^(2, R) 
with real eigenvalues, and does not use the fact that the monodromy matrices 
come from a specific second order self-adjoint eigenvalue problem. 

The next two results exploit the fact that the monodromy matrices are de- 
rived from Schrodinger operators to show positivity/monotonicity of various 
quantities defined above. The first type of result is a Sturm oscillation type 
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result for the energy variable E, and shows that a±,5±, (3± are of a fixed sign. 
Since, by the above calculation, these parameters represent a sense of rotation 
of the eigenvectors, this result shows that (as a function of the energy) one 
eigenvector winds about the origin in a positive sense as the energy increases, 
and the other winds in a negative sense. 

Proposition 1. The monodromy matrices M{E),N{E,x) satisfy the differen- 
tial equations (in E) 

Me ^ MJ^E) (19) 
Ne = Nje{E,x) (20) 

where the matrix $ is positive definite, and Q{E, 0) = and Q is positive definite 
for a; > 0. 

Proof. This calculation follows that given in Magnus and Winkler [22]. It is 
easy to see that Ne satisfies the differential equation 

N,E = JHNe + JHeN Ne{E, 0) = (21) 

and a solution by variation of parameters gives the formula 

Ne^ N [ N-^JHeN. (22) 
Jo 

Using the fact that N^JN — J the above becomes 

Ne = NJ [ N^HeN. (23) 
Jo 

It is clear that if He is positive semi-definite then / N^HeN is also positive 
semi-definite, and is in fact positive definite for x > Q unless there exists a non- 
zero vector V independent of y such that HeN{E, y)v = for all y. For the 
equations considered here it is easy to see that no such v can exist, and therefore 
N^HeN{E, y)dy is positive-definite. Similarly M{E) has the representation 

Me = MJ f U\E, y)HEU{E, y)dy. (24) 
Jo 

□ 

Remark 2. The fact that Me — MJ<^ (and similarly Ne) is simply the fact 
that the Lie algebra associated to the symplectic group is matrices of the form 
JH with H symmetric. The positivity of $ is essentially the Sturm oscillation 
theorem. 

Corollary 1. Let v±,w± be the eigenvectors in a gap for the periodic and 
defect monodromy matrices respectively, and a± ,I3±,S± be defined as in Eq. 
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dv± 



- a±Jv± (25) 

^ = P±J^± (26) 
dw± 
dx 



S±Jw± (27) 



then /3+,q;+ > and /3_,Q!_ < 0. Further, if the Hamiltonian matrix in the 
deject region Hj^f,f(E, x) is positive-definite, then (5+ > and S- < 0. 

Proof. We know that the /3± satisfy 



P+ = 



A+(w_, J Ne w-) 



1-Xl 

Using the fact that Ne = NJQ and N^JN = J we have that 

X+{w+,JNjew+) 



-X+{w+,N-*Qw+) 
{w+,Qw+) 



Ai-1 



> 



and similarly /?_, a+, a_. Using the fact that O is positive-definite and |A+p > 1 
we have that a+ > 0, < 0. 

The second part follows similarly: JNx = —HdefN, thus 



<5± = +^ +r 72 (28) 



1-A^ 

with Hdef positive definite, and the above argument follows 



□ 



3 Counting the Defects 

In this section we present the main results, which allow us to estimate the 
number of gap eigenvalues created by a defect potential in terms of the winding 
number of the monodromy matrix at the band edge. Again the basic idea is to 
make a homotopy argument in x of the generalized Evans function /(i?, x), and 
count the eigenvalues as the are created. 

The first important result is the following: 
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Lemma 2. Suppose that E G TZdef- The Evans function f{E,x) and its first 
partial with respect to the energy ^ cannot vanish simultaneously. For energies 
in the classically allowed region for the defect ^ cannot vanish at a zero of the 
Evans function, and has the same sign as ^ . 

Proof. Obviously we have the following expression for fs ■ 

fE = {%,JNv+} + {v-,JN^) + {v-,J% v+). (29) 

Using the identities ^ = A'' JO and N^JN = J the above expression reduces 
to 

fE = a^{v^Nv+) - a+{v^N-*v+) - (w_7V-*0w+). (30) 

We would like to compute the signs of the various terms above at a zero of the 
Evans function /. Note that, at a zero, we have Nv+ = iiv-, or N~^v- = 
lj,~^v+. Using this fact we have 

CK-i- 1 

fs = m- {v+,Qv+) (31) 

Since is positive definite, q;_ < and > we have that cannot vanish 
at a zero of the Evans function and has the opposite sign from (i. A similar but 
more straightforward calculation gives 

= {v_,JN,v+) = {v_,JJHa,sNv+) (32) 

= -fl{V-,HdefV-) (33) 

where the last step assumes that one is at a zero of f{E,x). Again iJrfe/ is 
positive definite in the Schrodinger case for energies in the classically allowed 
region (for the defect) and in the other cases without restriction. □ 

This result has a small technical drawback, in that the derivative of the 
Evans function with respect to energy diverges near the band edges of the pe- 
riodic problem, since the spectral quantities have square-root type singularities 
near the band edges. These can easily be eliminated by choosing appropriate 

coordinates in the gap. 

Lemma 3. For a given gap Gj = (-E'2j+i, £^2^+2) C Tiper define the new variable 
E by the invertible map 



E : 



dE 



ith iVj = Jj^ 

function of E in [0,0;^], and f and f^ cannot vanish simultaneously in [0,a;. 



that maps [i?2j+i, £^2^+2] to [0,0;^], with toj = J^^'^^ ^Jk^E) i ' f is aC^ 
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Proof. This is a straightforward chain-rule argument. The eigenvalue equation 
— k{E)\ + 1 = is equivalent to 

- 1 = ±\±y/k^{E)-4. 

From this together with the relation 

dE _ 1 

dE ~ yPpy34' 

we have that 

f^{E) = - -^l{^^'v^,'^>v^)^-(M-'v+,'Pv+)--{v+, ev+)JkiE)^ -4. 

dE M ^J■ ^ 

At the band edge, it holds k{E)'^ — 4 and AIv± — X±v± with |A±| ~ 1, and 
v+ (0) = V- (0) so it follows that 

lim f^iE) = -tl±l{v+{0),^0)v+iO)) 

where f+(0) is the (unique) eigenvector of M at the band-edge E' = 0, /i is 
defined by iV(0, x*)u_|_(0) — /iti+(0) which holds at a zero of /, and $ is the 
positive-definite matrix defined in (|20p . A similar result holds for the upper 
band-edge E = ujj □ 

Remark 3. The above is the hyperelliptic integral associated to the spectral 
prohlem \2Ji\ \3l . Note that the integral is necessarily convergent, since the Floquet 
discriminant k{E) necessarily crosses ±2 transversely at the edge of a non-empty 
gap. 

This lemma has the following obvious implication 

Corollary 2. Suppose that the gap (G) — (0,Wj) is a gap (parameterized by E 
as in lemma\^. A given zero level set of the generalized Evans function f{E, x) 
is given by a function E{x) which is either defined for all x e [0, 1], or leaves 
the range (0,0;^) at some point(s). In the classically allowed region E > q^ef 
this function is strictly monotone and thus invertible, x{E). 

Proof. A standard implicit function argument. Local existence is obvious. Since 
the Evans function and its derivative with respect to energy cannot vanish simul- 
taneously we have uniform control on the Lipschitz constant on [0,(jjj] x [0, 1], 
and thus the level set can be constructed for all x as long as E{x) remains in 

□ 

From this it should be clear that in the classically allowed region for the 
defect the number of roots of the Evans function f{E, 1) can be expressed in 
terms of the number of roots of the Evans function at the band edges and 1 . 
This is the content of the next proposition. 



11 



Proposition 2. Suppose that the gap {Q^ujj) is in the classically allowed region 
for the defect potential. Let N be the number of roots of f{E,l) in the gap 
E s (0, 1) rii be the number of roots of f{0,x) for x € (0, 1] and n2 be the 
number of roots of f{u)j, x) for x G (0, 1). Then N = n2 — ni + 1. 

Proof. This is a fairly standard lemma - that simple roots persist under a ho- 
motopy that fixes the endpoints. 

Consider the region {E, x) G [0, uij] x [0, 1]. It is easy to see that for a; = the 
Evans function f{E, 0) is nonzero for E G (0, cjj) and vanishes at the boundaries 
E = {O^Uj}, since the monodromy map N{E,x) is the identity for x = and 
v± are non-degenerate on the interior of the band and degenerate at the band 
edges. Since x{E) is decreasing the zero of f{E, x) at {E = 0,x = 0) does not 
extend into the interior. If f(E.x) vanishes at the band edge (wj, 1) this zero 
also does not extend into the interior. By the above lemmas all other zeroes of 
the generalized Evans function can be continued to a zero level-set x{E) that 
is a monotone decreasing function. Thus, each zero of the generalized Evans 
function on the right-hand boundary E = ujj extends to a unique curve that 
must intersect the left-hand boundary = or the top boundary x =1. Thus 
we have n2 -|- 1 = A'' -|- m . 

□ 

The first homotopy argument shows that the number of roots of the Evans 
function (in E) is equal to the difference in the number of roots (in x) on the 
band edges of the generalized Evans function. We would next like to be; able 
to estimate the number of roots of the generalized Evans function in terms of 
the band-gap structure of the defect fundamental solution matrix. The next 
lemma shows that zeros of the Evans function on the boundary can only occur 
in gaps or band-edges for the fundamental matrix N{E, x) for the defect mode, 
but never in bands. 

Lemma 4. Suppose that the generalized Evans function f{E,x) vanishes for 
an energy at a band-edge of the essential spectrum, {E,x) = (0,a;*) or {E,x) = 
iLOj,x*). Then kl^f{E,x*) > 4. 

Proof, li E = {] OT E = Ldj the eigenvectors w+ and are linearly dependent, 
and thus the condition that the Evans function implies that 

f{E, x)=(i^{v+, JNv+) =0^ Nv+ = I1V+. 

Thus Vj^ is an eigenvector of N{E,x*). Since v+ is a real vector it can be 
an eigenvector of A'^ only if A'^ has real eigenvalues, which is equivalent to 
klfiE,x*)>4. □ 

The above shows that zeros of the Evans function on the boundary can only 
occur in gaps of N. We would like to show that exactly one zero of the Evans 
function emerges from each gap. The next lemma, which will be needed to 
show this, says that gaps in x behave qualitatively like gaps in E for energies 
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in the classically allowed region, in the sense that in an open gap the Floquet 
discriminant necessarily crosses ±2 transversely. The analogous statement for 
gaps in E is standard (see Magnus and Winkler) and will not be proven here. 

Lemma 5. Suppose that xq is a periodic or anti-periodic point for the de- 
fect, kdef{E,xo) = ±2 and that the energy E is in the classically allowed re- 
gion for the defect. Then either the Floquet discriminant crosses transversely, 
^^^(-E, Xq) ^ 0, or xo is a double point N{E, xq) = ±/ and kdef = ±(2 — c(a; — 
xo))^ + o{{x — xof)) with c> 0. In particular for an open gap the eigenvectors 
are degenerate only at the band edges. 

Proof. We prove for the case kdefiE^xo) = 2. The case kdef{E,Xf)) ~ —2 
follows similarly. If kdef{E,XQ) ~ 2 then N takes the Jordan normal form 
N = I + K V (E) Jv where v is the (non-zero) eigenvector of N. Note that 
K can vanish, in which case N is the identity matrix. Taking the identity 
tr(iV2) - (tr(iV))2 = -2det(iV) = -2, differentiating and substituting into the 
Jordan normal form gives = «; < vJN^v >= —k < vH^efV > . Since 

the energy is in the classically allowed region for the defect the Hamiltonian 
is positive definite and the vanishing of t'r{N^) implies that k = and thus 

that the Jordan normal form is the identity. In this case we find that ^^i' = 
tr{JH{E, xq)JH{E, xq)) = -2 det{H{E, xq)) < 0. □ 

The next lemma shows that the Evans function has exactly one zero in a 
gap (in x) of N{E, x) if the Hamiltonian ifde/ is positive-definite. 

Lemma 6. Suppose E is a fixed energy such that Hdef{E, .t) is positive-definite, 
v+ a fixed real vector, and {xiow, Xhigh) is an interval such that 

• ^de/(^. Xlow) = 4 = kj^fiE, Xhigh) 

• ^defi^^ X) > 4:forX G {X low, Xhigh) ■ 

In other words {xiow,Xhigh) is a gap in x of N{E,x). Then the condition 

{v+JN{E,x)v+} = 
has exactly one root in x G [xiow,Xhigh]- 

Proof. This follows from a simple monotonicity argument. We assume that 
the gap is not a double point, in which case the lemma is trivially true. The 
eigenvectors of N{E, x) are degenerate at the band edges and (by the previous 
lemma) non-degenerate in the band interior. We can choose the eigenvectors 
w± in such a way that = W- at the lower band edge xiow From Corollary 
1 it follows that ui+ rotates clockwise and W- rotates counterclockwise until 
W-i^ = —W-, which (by the previous lemma) occurs at the upper band edge. By 
continuity it follows that the angle between w+ and W- goes through tt radians. 
The eigenvalue condition is equivalent to either w+ = v+ or w- = v+.lt follows 
that either v+ is equal to a band-edge eigenvector of N{E, x), in which case the 
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Figure 1: The picture of the band-gaps in the {E,x) plane. The bands can 
touch at a double point, but cannot cross. Gaps, however, can disappear and 
reappear. 



above eigenvalue condition has a simple root at the appropriate band edge, or 
is not a band edge eigenvector, in which case there is exactly one root of 
either w-|_ — or w- = v+ (but not both). This is illustrated in Fig. (1). □ 

Lemma |4] showed that the eigenvalue condition on the boundary of a gap in 
the essential spectrum could only have a roots in the set k'^^j:{E,x) > 4. The 
above lemma shows that the eigenvalue equation has exactly one root in each 
interval of this form. 

At this point we have shown that, for energies in the classically allowed 
region for the defect, eigenvalues persist under homotopy in x. We have also 
shown that there are certain distinguished intervals in x, [xiowt ^high]i where 
eigenvalues can enter the gap through the band edge E = E2j+2 or leave the 
gap through the band edge through the band edge E = i?2j+i- Further we 
have exactly one eigenvalue entering or leaving the gap in each such interval. 
Thus the number of eigenvalues in the gap is equal to number which enter 
through E = £'2^+2 minus the number which leave through E2j+i. This is not a 
particularly convenient charactertization, since it requires us to count intervals 
in the homotopy parameter x. Thus we give a second homotopy argument to 
count the number of gaps (in x) of N{E,x) in terms of the number of gaps (in 
E) for N{E, 1). 

Lemma 7. For 6 sufficiently small the band gap structure of N(E, S) in E for 
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E in a fixed interval about q{0) consists of 

• Gap for E < q{0) - o{l) 

• Band for E > q{0) + o{l) 

Further there is a unique E in {q{0) — o(l),g(0) +o(l)) where ty:{N{E,6) = 2, 

and the crossing is transverse. 

Proof. This result follows from Taylor expansion, self adjointness, and some ele- 
mentary analytical considerations. First note that N{E, 0) = /. Differentiating 
gives tvN,:,,iE,0) = -2{E-q{0)), trN^EiE,0) = and trNEEiE,0) = 0. So 
for smalla; A;(i;,a:) = triV(£^,a:;) w 2- {E-q{0))x^ + 0{x^). Therefore, it is clear 
that for X sufficiently small, k{E, x) > 2 fov E < q{0) - 0{6) and k{E, x) < 2 
for E > q{0) + 0{S). Next we need to show that there exists a single transverse 
crossing in (g(0) - 0{6), g(0) + 0{S)). To do this note that N{E, x) is the mon- 
odromy map for a self-adjoint operator, and standard results (sec Magnus and 
Winkler) show that the only possible band edges are a proper band edge, which 
crosses transversely (in this case k{E,x) = 2,kE{E,x) < 0) or a double point 
(in this case k{E, x) = 2, ksiE, x) = 0, kEsiE, x) < 0). Thus the first crossing 
must be transverse. To show that this crossing is unique for sufficiently small 
S we note that self- adjointness implies that the Floquet discriminant cannot 
have a critical point inside of a band. Suppose that there were two crossings in 
{q{0) — 0{5), g(0) + 0{S)). Then, since critical points inside of a band are disal- 
lowed, there must be a gap in {q(0)—O{S), q{0)+O{5)), and thus two antiperiodic 
points. Thus the Floquet discriminant must go from +2 to —2 in a distance of 
0{S), and we have an estimate of the form \k{Ei, S) — k(Ei + 6, S)\ > 0{S^^). 
Since the Floquet discriminant is continuous this cannot hold for 5 arbitrarily 
small, and thus the crossing is unique for small S. □ 

Now we arc in a position to prove the next result: 

Proposition 3. Suppose that E is in the classically allowed region for the defect 
mode. The number of bands of N{E, x) in x G [0, 1] is the equal to the number of 
bands of N{E, 1) for E <E (— oo, i?] . The number of gaps of N{E, x) in x & [0, 1] 
is one less than the number of gaps of N{E, 1) for E G {—oo,E]. 

Proof. This again makes use of a homotopy argument. The basic observation 
is that, while the gaps can collapse to a point and thus cannot be continued by 
an implicit function argument the bands can be so continued. 

It is somewhat more convenient to work on the the interval [6, 1] x [e, E] to 
avoid trouble with noncompactness in E and the degeneracy at a; = 0. Since 
qdef is assumed to be piecewise we choose e < infxeiQ.^li^)- fhis case it 
is clear from the obvious lower bounds that there are no bands for energies less 
than or equal to e. From the above lemma we can choose 6 sufficiently small 
such that E < q(Q) — o(l) is a gap, E > q{Q) + o(l) is a band, and there is a 
single periodic point a,t x = S,E = g(0) +0(1). Since the bands are by definition 
open we can chose 6 such that the number of bands for x G [0, 1] is the same as 
the number of bands for x e [5, 1]. 
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We define the bottom boundary B to be [e, E] x S, the top boundary T to 
be [e, -E) X 1, the right boundary i? to be ii^ x [6, 1], and the left boundary B to 
be e X [S, 1] . Since ^ cannot vanish on the interior of a band it follows from the 
implicit function theorem that a point in a band defines a local level set of the 
the Floquet discriminant E{x). Since we have global control of the Lipschitz 
constant this level set is defined for all x. 

By construction the left boimdary L is chosen to be below the essential 
spectrum and has no bands. The bottom boundary B has a single band (g(0) + 
0{5), E), and by construction dk/dx < on this band. This band continues onto 
the right boundary R. Since E is in the classically allowed region dk/dx < 
for all bands on top boundary R. Thus all bands of the left and top boundaries 
move into the interior transversely as x increases. Since a level set of k{E, x) on 
the interior of a band is defined globally as a function of x (again denoted E{x)) 
each point on a band interior on the top boundary moves into the interior as 
X is increased, and can be followed until it intersects the right boundary. This 
defines a map of bands on the right boundary to bands on the top boundary. 
The converse argument defines the inverse map of bands on the top boundary 
to bands on the right boundary. Since we have an invertible map of band on 
the top boundary to bands on the right boundary the number of bands on the 
top boundary and on the right boundary must clearly be the same. Since the 
bands are open and disjoint the complements must also have the same number 
of components. The bottom boundary contains one gap, the semi-infinite one, 
and the remainder of the gaps occur on the right boundary, whence the gap 
count. 

□ 

Remark 4. Note that the fact that the derivative dk/dx is of the same sign 
as dk/dE is critical to this argument, and depends crucially on being in the 
classically allowed region for the defect. Without control on the sign of dk/dx (or 

equivalently dx/dE) there is no guarantee that a band could not enter and leave 
the region repeatedly. Physically this seems to be a non-resonance condition. 

From this it follows that we have the following bounds on the number of 
defect eigenvalues. 

Theorem 1. Suppose that a gap Q is in the classically allowed region for the 
defect mode. Define the set 

nd,f = {E\{tv{M{E))f>A} 

and the integers ng and ngg to be the number of connected components of the 
sets TZdef n Q and TZdef H dQ respectively. Then, the number of defect modes 
Dg, satisfies the inequality 

ng + 1- noG <Dg<ng + l. (34) 

Proof. We'll prove the case ngg = first, in which case the count is exact. If 
noG = 0, i.e., neither edge lies in a gap for the defect monodromy matrix, then 
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from the previous proposition and lemma 8 it follows that the number of gaps of 
N{E, x) in X is the equal to the number of gaps of N{0, x) in x ioi x G [0, 1] plus 
no- But, since in each gap (in x) there is a unique eigenvalue, we obtain that 
(in the terminology of Proposition 1) n2 — ni = no, and therefore Dq = riQ + 1. 

On the other hand, for each of the edges 0, 1 that lie in a gap for the defect 
monodromy, the defect modes can be reduced by 1. This is true, since the zero 
level curves for the Evans function might not exit the x— interval [0, 1) in the 
case where E'jotu lies in a gap, or the defect mode might be generated outside of 
[0, 1) in the case where Ehigh hes in a gap. □ 

Remark 5. It is worth making a number of remarks about this result. First we 
note that one could easily find upper and lower bounds on the number of defect 
eigenvalues generated in a gap by counting Dirichlet eigenvalues of the defect. 
However using this method there seems to be no clear criterion for guaranteeing 
that the count is exact, which will be necessary for the main result of this paper, 
the large gap number estimate. 

As noted earlier the set Rdef is the union of the resolvent set and the periodic 
and anti-periodic eigenvalues of the periodic Schrdinger operator with potential 
given by the periodization of the defect potential. Thus the above result, roughly 
speaking, gives an estimate of the number of defect eigenvalues in terms of the 
number of connected components of the intersection of the resolvent sets. 

The above calculation can also be interpreted as a Maslov index calculation. 
Given a curve N[E) in the symplectic group and a Lagrangian subspace C (in 
this case any one dimensional subspace ) the Maslov index is a signed count of the 
number of intersections of N{E)C with C In our case lemma ?? guarantees that 
the index is always positive, so the Maslov index actually counts intersections. 
This count depends, of course, on the Lagrangian subspace chosen but in the 
special case ngg = the Maslov index fic^s the same for all Lagrangian subspaces 
and the number of defect eigenvalues is one plus the Maslov index of N{E) for 

E £ {E2j+l,E2j+2)- 

Finally we remark that, in the case of defect potentials which are constant 
the set TZdef consists of single points. In this case it is easy to see that the count 
is exactly ng + 1 regardless of whether or not ngg = 0. 

3.1 Large gap number asymptotics 

In this section we present results for the high gap-number behavior of the number 
of defect eigenvalues. It is known [7j that the n-th Dirichlet eigenvalue for the 
problem 

—Uxx + V{x)u — A„u, m(0) — u{L) — 
behaves asymptotically as 

n^TT^ 1 

A„ + T V{x)dx + o{l). 

^ ^ Jo 

We have established previously bounds for the number of eigenvalues in each 
gap. In particular we have that if the gaps of the (periodically extended) defect 
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problem do not intersect a gap of the periodic problem, then that gap has exactly 
one defect cigcnvahic. Further we know the following standard facts: 

• The n*^ Dirichlet eigenvalue is contained in the n*'' gap. 

• The width of the n*'* gap goes to zero for large n. In particular iiV £ L2 
then the sequence of widths is in I2, with stronger decay estimates if V 
has additional smoothness properties. 

Thus, if we denote the n-th Dirichlet eigenvalue of the periodic problem by /i„ 
and the m-th Dirichlet eigenvalue for the defect problem by /i„ then we can 
expect to have one eigenvalue in each gap, unless Hn ~ Mm) for some integers 
m, n. To quantify this we first define the following set of exceptional energies, 
which detects possible overlaps between the Dirichlet spectra of the periodic 
potential and the defect potential: 

Definition 1. For a real number a we define the set Ta as follows: a number 
y belongs to Ta if for every 6 > there exists an infinite, strictly increasing 
sequence of pairs of integers {Nk,Mk}'^i such that 

\y - Mk{Nk - Mka)\ < 5 

In other words, the set !Fa represents the (appropriately scaled) asymptotic 
distribution of errors in rational approximations of a. In the next lemma we 
note some simple properties of the set J^a- 

Lemma 8. The set Ta has the following properties: 

• The set Ta is never empty. If a is rational then Ta contains only the point 

y = o. 

• If a is irrational then there are at least countably many points in Ta ■ 

• The set Ta is invariant under the modular group PSL(2, Z); if y € Ta 
with Mk{Nk — aMk) — > y and a = ^^^^^^ with n\ni — 112123 = 1 then 
yeFb. 

Proof. The first property is obvious, since it is easy to see that a given non-zero 
integer can be written as a difference of squares in at most a finite number of 
ways. For the irrational case note that from elementary number theory if a is 
irrational there exists an infinite sequence of pairs {N, M) such that — a| < 
jj^j-jthus an infinite sequence of pairs {N,M) such that M{N — Ma) € (—1, 1), 
and by compactness a limit point. Note that if y G Ta then j^y G Ta for all 
integers j, so the number of points in Ta is at least countably infinite. 

Finally to see the invariance of Ta under the modular group note that a 
straightforward calculation shows that if the sequence {N^, M^) satisfies Mk{Nk— 
aMk) y then the sequence Nj^ =niNk — n2Mfe, M'y. = n\M}. — nsNk satisfies 
KiK - 6M^) ^ y □ 

Finally, we characterize the set Ta for some classes of real numbers a. 
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Proposition 4. Suppose a has the continued fraction expansion a= [ao,ai,a2 ■ . ■ 

• If {aj}'j°^Q is eventually periodic (in other words if a is a quadratic irra- 
tional) then Ta is a discrete set of point. 

• If {a2j}j^Q is unbounded then Ta contains the negative half-line (—00,0]. 

• If {a2j+i}jtQ is unbounded then Ta contains the positive half-line [0, 00). 

• V WA'jLi bounded then Ta does not contain some interval about the 

origin. 

In particular for Lesbesgue almost every a both the even and the odd terms in 
the continued fraction expansion are unbounded and thus we have Ta = R, the 
whole real line. 

Proof. In the case of quadratic irrationals it is relatively easy to compute ex- 
plicitly what Ta is. Let f{x) = rixx^ + n2X + ns with ni,2,3 S ^ relatively prime 
be the quadratic polynomial with root a. Then Ta consists of all numbers of 
the form = . ^ , where i is any integer which can be represented in 

J y"'! y/n^— 4,711713 

the form niN"^ + n2NM + n^M'^ = j. 

To see this note that standard results in the theory of quadratic Diophantine 
equations show that the existence of one solution guarantees the existence of a 
family of solutions. Let Nk,Mk be increasing sequences of integers such that 
Nk/Mk — » a. Obviously niN^ + n2NkMk + n^M^ = j for some integer j. Since 
/(a) = 0, it follows by the mean value theorem that 

/(a) - f (Nk/Mk) = f{xo)ia - Nk/Mk), 

for some Xq between a and N^/Mk, and so we obtain f'{xo){a — Nk/Mk) = 
-3 /Ml. Therefore 

Mk{Nk - aMk) = 

f [xo) 

for xo & (a, Nk/Mk). Thus all elements of Ta are of the stated form. To see that 
all such numbers are actually arise note that the existence of one solution to 
niN"^ + n2NM + n^M^ = j implies the existence of a family of such solutions, 
which by the above must satisfy 

Mk{Nk - aMk) = 

f [xo) 

To see the second and third claims we consider the following doubly indexed 
sequence: Nj^k = jPk,Mj^k = jQk, where ^ arc the continued fraction ap- 
proximants to a. Standard results show that |P| — a^Qfel ^ so if Ok is 
unbounded there exists a subsequence such that — Q^Q^ = ^ 0. Further 
Ck is negative (resp. positive) if k is even (resp. odd). Note that for a fixed 
integer j we have that j^Pk — oPj'^Qk = j'^^k- It is clear that for any y having 

the same sign as Ck if we take jk = [^J we have | j^e/j — y\ < )• 
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The fourth observation foUows from the weU-known fact that a number with 
a bounded continued fraction cannot be approximated by rationals to better 
than quadratic order: if the continued fraction coefficients satisfy Ui < m then 
one has a lower bound of the form 

\N/M-a\ > f{m)/M'^ 

with /(to) > 0, from which it follows that the interval \y\ < f{m) is not in the 
set J^a- 

The last assertion follows from a trivial modification of the proof in Hardy 
and Wright that the set of numbers with bounded continued fraction coefficients 
has measure zero. □ 

Remark 6. The case which is still incompletely understood is that for which 
the continued fraction expansion has bounded coefficients but is not periodic. It 
would be interesting, although likely very difficult, to classify JF^ based on the 
distribution of the continued fraction coefficients. 

Now we are in a position to state the main result, which relates the asymp- 
totic number of eigenvalues in a gap to the properties of the set J^a- 

Theorem 2. Assume that the widih of the defect is norm.a,lized to 1, and, as be- 
fore denote the period of the periodic potential by a. Define the energy difference 
Aq to be the mean of periodic potential minus the mean of the defect potential: 



Ag= / qdef{y)dy-- [ qper{y)dy. 
Jo 0. Jq 



Then: 



• // the energy difference does not belong to the set of exceptional energies, 

^ J-'a, then every sufficiently large numbered gap contains exactly one 
defect eigenvalue. 

• // e then there exists a sequence of exceptional ga,ps. Every suf- 
ficiently large numbered gap which is not in the sequence of exceptional 
gaps contains exactly one eigenvalue. Every sufficiently large numbered 
gap which is in the sequence of exceptional gaps contains 0, 1, or 2 defect 

eigenvalues. 

Proof. Denote the n-th Dirichlet eigenvalue of the periodic problem by /Lt„ and 
the TO-th Dirichlet eigenvalue for the defect problem by fim- Then, 

fJ-n-Jj'm ^2 ■^^'^^ - + 0(1) (35) 

'-^ ^ (n + am)(n — am) — Ag + o(l). (36) 
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For n, m large this is clearly bounded away from zero unless n = am + o(l), in 
which case this becomes 



27r2 

— Mm ~ "'^^ — Ag + 0(1). (37) 

If Ag does not belong to the set of exceptional energies J-a then this quantity is 
eventually uniformly bounded away from zero. Thus all sufficiently large Dirich- 
let eigenvalues of the periodic problem and the defect problem are bounded away 
from each other. Since the Dirichlet eigenvalues are contained in the gaps, and 
the width of the gaps is approaching zero, it follows that the gaps are eventually 
non-intersecting. From the count in Theorem[T]it follows that each gap contains 
exactly one eigenvalue. 

If AE £ Ta then it is clear that at most one of the gaps of the defect problem 
can intersect a gap of the periodic problem. Hence the gap of contains either 
0, 1 or 2 defect eigenvalues according to Theorem [TJ 

□ 

To close we construct an example to show that the the exceptions considered 
in the previous theorem are real. In particular we construct a combination of 
periodic and defect potentials such that a particular (infinite) subsequence of 
gaps contains two defect eigenvalues. 

Example 1. In order to construct this example we would like to choose the 
periodic and defect potentials in such a way that the gaps in the essential spec- 
trum of the periodic problem are comparatively wide. For this reason we choose 
a piecewise constant (Kronig- Penney) potential. We choose the period to be a 
quadratic irrational since this guarantees that we have a clean description of the 
set Ta- For reasons of tradition we take a = 4> the golden mean, although any 
quadratic irrational would serve. We take the defect potential to be constant, 
Qdefix) = Qdef, SO the Set TZdef consists of a union of points, and the eigenvalue 
count is guaranteed to be exact. In particular for the periodic potential we take 
the Kronig-Penney potential 

r -A for 0<x< 0/2 
qper(X)^<^ ^ /or0/2<x<0 

where A is a constant to be determined later. The gaps grow in size with A, 
and our strategy is to choose A sufficiently large so that appropriate points of 
TZdef cLre eventually contained in the appropriate gaps, as this will guarantee 
two eigenvalues. 

The transfer matrix for ~Uxx + qper{x)u — Eu can be easily calculated to be 
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From this, the location of the band- edges is given by roots of the equation 



-{^ + m) --(^) X sin(^) . ±2 . 

From some straightforward asymptotic analysis it can be seen that the loca- 
tion of the j*'' gap is approximately 



OU-'l ^ + ^ + Oir')) for J » 1 odd 



4^ + 0(j-^), 2^ + ^ + 0{r')) for J » 1 

Since the defect potential is constant the periodization is constant. Thus the 
gaps close to double points and the set TZdef is given by 



The set consists of real numbers of the form 



= j e ±{1,4,5,9,11,16,19,20,25,29,31,36...} (39) 

where the sequence is all positive integers j which are representable in the form 
j = — nm — (sequence A031363 in Sloane's encyclopedialSSi). The mean 
of the periodic potential is zero, so if the mean of the defect potential is chosen 
to be in the set ^^J'cj) then there is a infinite sequence of gaps which potentially 
have more than one defect eigenvalue. We will, somewhat arbitrarily, focus on 
the j = 11 term of T^. This term is interesting since it admits two distinct 
families of solutions to the Diophantine equation nf, — n^mfc — — 11. These 
two families are given by the following sequences: 

mfe = {1,5,14,37,97...} = {4, 9, 23, 60, . . .} (40) 

mfc = {2,7,19,50,131...} = {5, 12, 31, 81, . . .}. (41) 

Each of the sequences nk,mk satisfies the recurrence ak+i = 3ak — flfc-i and 
i^k ~ fikTUk — m\ — 11. It is straightforward to compute that the nk,mk satisfy 
the following asymptotic relation 

1 + V5 11 121 , _4^ 

rukiuk — TOfc) = — = , ^— + 0(mj^ ). 

^ V5 m^vl25 

(This is easiest to see if one notes that — Ukmk — rnf. — 11 can be factored 
over Q(\/5) as {nk — (t)mk){nk + 4'~^rnk) — 11 and proceed from there.) 

We look at the Gn^ , the n^^ gap in the essential spectrum with nk an element 
of one of (|40|) or ()4ip . and the m}j^ point in TZdef , which is located at Tr^m^ + 
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qdef- The distance between the center-point of ^.n-d T^'^'n^k ~^ Idef is given by 

2 2 , 2( , "-feN/ "-feN 

nm^. + qdef = t: [mk + —j(mk -) 

0^ (p (p 

227r2 12l7r2 _2, 

~ (Idef J=- + 7-5 2 + "('^/c )■ 

Choosing qdef = (again this is the j = 11 element of J- c/,) causes the 

leading order terms to cancel. Given the gap asymptotics in i38\) it is clear that 
for Uk odd the point /im^ is always eventually contained in the gap Qn^ , since 
the width of the gaps Qn^ decays more slowly that the above error. For mk even 
the point fimk will eventually be contained in the gap Qn^ */ ^ chosen to be 
sufficiently large. If we choose 

3cj}^A^ 121^2 



> 
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or equivalently 

then the radius of the Uk gap is asymptotically larger than the distance between 
the ruk point in TZdef cind the center of the gap, and thus the ruk point in TZdef 
is contained in the nj!^ gap. Thus ng — l,nQg = and the gap Qn^ contains 
precisely 2 eigenvalues. 

Note that this example is fairly robust. For any positive A the odd gaps Qn^ 
for Hk in the sequences (|40p or (|41|) get exactly two eigenvalues: the lower bound 
on A is necessary only to insure that the even gaps (which are asymptotically 
narrower) get two eigenvalues. Also note that the choice of qdef (x) is not par- 
ticularly important: as long as qdef is smooth enough that the intervals in TZdef 
decay at least as fast as the gaps in the essential spectrum then (for sufficiently 
large contrast A) the above example can be made to work. Similarly, as noted 
above, the above construction goes through for any quadratic irrational ( in fact 
<f) is in some sense the worst example since it is the least well approximated by 
rationals). For a general irrational a and general potentials qdef{x) and qper{x) 
there is a somewhat delicate interplay between the size of the intervals in TZdef 
and TZper and the error in the rational approximations of a, and while it is clear 
that the sequence of extraordinary gaps can contain at most two eigenvalues it 
is difficult to determine the exact number. 



4 Conclusion 

In this paper, we analyzed in detail the emergence and distribution of defect 
eigenvalues in the gaps of the essential spectrum of the linear one-dimensional 
Schrodinger equation, with a potential consisting of a periodic part plus a com- 
pactly supported defect part. We used an Evans function technique that reduces 
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the problem to that of finding the zeros of an analytic function, and by means 
of a homotopy argument we were able to count the eigenvalues as they emerge 
from the band edges. It is found that if a gap is in the classically allowed region 
for the defect potential, then the number of defect modes in the gap can be ex- 
pressed in terms of the number of roots of the Evans function at its band edges. 
It is shown that for a gap in the classically allowed region for the defect poten- 
tial, the number of zeros of the generalized Evans function is no larger than the 
number of gaps of the defect problem in terms of the homotopy parameter x. 
In addition, this number was found to be the same as that of the gaps in terms 
of the eigenvalue parameter E. As a result, bounds for the defect modes are 
given in terms of the gaps of the defect monodromy matrix that intersect the 
gap. Finally, we proved the following significant generalization of Zheludev's 
theorem: the number of point eigenvalues in a gap in the essential spectrum 
is exactly 1 for sufficiently large gap number unless a certain Diophantine ap- 
proximation problem has solutions, in which case there exists a subsequence of 
gaps containing 0, 1, or 2 eigenvalues. We stated some conditions under which 
the solvability of the Diophantine approximation problem can be established, 
and we included an example where a particular (infinite) subsequence of gaps 
contains two defect eigenvalues. 

There are a number of interesting open questions. It would be interesting to 
understand the structure of the set J^a for a irrational with bounded coefficients 
but not a quadratic irrational, although this would likely depend in a sensitive 
way on the distribution of the continued fraction coefficients. Also, since the 
bounds are expressed in terms of the intersection of resolvent sets it would be 
interesting to understand how these are effected by isospectral flows on the 
potential(s), which obviously leave these sets invariant. We believe, although 
we have not yet checked this, that by flowing the periodic potential according 
to the Korteweg-DeVries hierarchy (which obviously leaves TZper invariant) one 
should be able to force any particular gap to achieve any of the possibilities 
allowed by ([M]). 
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